RADIAL STANDING AND SELF-SIMILAR WAVES FOR THE 
HYPERBOLIC CUBIC NLS IN 2D 



PANAYOTIS KEVREKIDIS\ ANDREA R. NAHMOD^ AND CHONGCHUN ZENG^ 

Abstract. In this note we propose a new set of coordinates to study the hyperbolic 
, or non-eUiptic cubic nonhnear Schrodinger equation in two dimensions. Based on 

^SJ ' these coordinates, we study the existence of bounded and continuous hyperbohcally 

radial standing waves, as well as hyperbohcally radial self-similar solutions. Many 
' of the arguments can easily be adapted to more general nonlinearities. 



1. Introduction 

' The hyperbolic cubic NLS equation has been studied in [181 US 120] (see also [281 

. [121 [TT] and references therein) and has received attention by both the physics and 

the applied mathematics communities. This model has been used for examining the 
evolution of optical pulses in normally dispersive (quasi-discrete) optical waveguide 
array structures [161l24j. as well as more generally in normally dispersive optical media 
^ |131ll4j. These studies have, in turn, motivated a number of works examining the type 

. of coherent structures, such as the experimentally tractable X- waves [IflHHlH], as 

I well as more elaborate structures such as dark-bright [21] or vortex-bright solitary 

I waves [17J. On the one hand, the spontaneous emergence of nonlinear structures such 

■ as the X-waves has been revealed by means of mechanisms such as the modulational 

. instability \9\ and their nonlinear dynamics has been considered jlO], but also more 

recently studies have begun to systematically address their spatio-temporal interactions 

m\. 



The equation arises when to the normally-dispersive (l-l-l)-NLS, effects due to the 
^ ■ diffraction (in the quasi-discrete setting) or anomalous dispersion (in the continuum 

. setting) are incorporated as an additional term. The (2-|- l)-dimensional quasi-discrete 

equation [HI [2l] is described as 



(1.1) - fS^ + ^^^"+1 + ^"-i) + = 

where U is the complex amplitude that represents an envelope of the relevant wave 
field, z is the longitudinal direction along which propagation occurs; and n and t are 
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the two transverse axes (which represent space and time respectively). In that same 
vein, it is relevant to mention that very recently also (3+l)-dimensional generalizations 
of the above quasi-discrete setting have emerged, giving rise to the so-called discrete- 
continuous X- waves in photonic lattices |22j . 

The continuous version of the (2 -|- l)-dimensional equation can be written as 

iut + Uxx — Uyy + 7|n|^n = (3HNLS) 

and is related to both the Davey-Stewartson and Ishimori systems \19\ 128], In fact, in 
[18J Ghidaglia and Saut showed that (3HNLS) satisfies the same linear Strichartz esti- 
mates as its elliptic counterpart and relying on them, they proved the hyperbolic-elliptic 
Davey-Stewartson is locally well posed in L^(R^) with time of existence depending on 
the profile of the initial data and globally well posed for sufficiently small data. The 
same argument proves that (3HNLS) enjoys the same well posedness results. In |20] . 
Ghidaglia and Saut showed that there are no nontrivial localized standing wave solu- 
tions to (3HNLS); more precisely they proved that there are no solutions of the form 
e'^'^^(j){x) with (f) £ n Hf^^ and any ^ G M (regardless of dimension). From the phys- 
ical point of view (3HNLS) or equation (|2.1|) below represents the canonical model for 
the study (experimentally and also numerically) of X- waves, as considered in [T3l [H]: 
see also the earlier and more recent reviews of \12\ lllj. It should also be mentioned 
that the linear counterpart of the equation has also been considered extensively [11^ [8] 
and has given rise to interesting recent experimentally examined structures such as the 
linear alternative to light bullets proposed in [7] in the form of Airy-Bessel wavepackets. 

The major differences between the regular Laplacian A = dxx + dyy and the hy- 
perbolic Laplacian Oxy = dxx — dyy include that a) the latter has nontrivial char- 
acteristics; b) the former is invariant under the Euclidean rotations on M? while the 
latter is invariant under the hyperbolic rotations; and c) the former corresponds to the 
positive Dirichlet energy and the latter to the energy E{u) := I'UxP — \uy\'^dxdy 
which is indefinite. Due to observation a), this PDE system may have singularities 
forming along the characteristics and thus we first find some explicit compatibility 
conditions along the discontinuities of weak solutions which are piecewise smooth. Ob- 
servation b) naturally motivates us to consider the hyperbolic coordinates (a, a), given 
by X = zbacosha and y = itasinha. Since the hyperbolic NLS is invariant under the 
hyperbolic rotations (a, a) — )• (a, a + ao), it is natural to start our study with solutions 
with the corresponding symmetry - continuous and piecewise smooth hyperbolically 
radial (or equivariant) standing waves and hyperbolically radial (or equivariant) self- 
similar wave£| - as our initial attempts to further our understanding of this PDE. 
They also usually play important roles in studying the dynamics of (3HNLS). In par- 
ticular, states which are hyperbolically radial share some similarity with the profiles 
of the two-dimensional analogs of the so-called X- waves observed in numerics; as a 
relevant example see e.g. Figs. 18.6-18.7 of [12| (cf. also |llj). Since the measure 
dxdy = \a\dada, non-zero hyperbolic radial states do not belong to L^(IR^). We do 



By self-similar waves we mean standing waves in the self-similar variables; c.f. (|3.6p (|3.7|) . 
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not view this as an extremely undesirable situation. On the one hand, some important 
elementary waves, like plane waves and kinks, do have infinite norm. On the other 
hand, since the hyperbolic Laplacian U^y does not correspond to a positive energy 
and has characteristics, it is reasonable, at least as an initial step, to relax the finite 
energy requirement and look for hyperbolically radial standing waves which may have 
asymptotic values other than zero at spatial infinity in each hyperbolic quadrant. 

In constructing these hyperbolically radial (or equivariant) waves, we encounter an 
ODE with a regular singular point. Though for this specific equation there are many 
standard results concerning the existence of local solutions, we provide a different proof 
based on the invariant manifold theory in dynamical systems which is a generalization 
of a trick used in [15]. Unlike the power series method or variational method, this 
approach requires minimal smoothness and no extra structure. Our aim in the present 
paper is to use such techniques and, especially, the hyperbolic coordinate system in or- 
der to construct, for the first time to our knowledge, nontrivial standing wave solutions 
to the (3HNLS) equation. In addition, we also consider radial self-similar solutions in 
the same coordinate system. The introduction of the hyperbolic coordinates and the 
analysis on hyperbolically radial waves represent the first step in a longer term program 
to further our understanding of PDEs involving the hyperbolic Laplacian operator. 

In the rest of the paper radial always means hyperbolically radial. 



In the following, we give the compatibility conditions satisfied by weak solutions 
following the standard argument. Let us denote by Uxy '■= d"^ — dy and for (x, y) G 
suppose u{t,x,y) is a solution of (3HNLS) with 7 = 1; ie. u satisfies: 



in the distribution sense and is smooth except on a surface T C M^^^, let us denote 
by Q± the two domains separated by T and by = {Nq, Ni, N2) the unit outward 
normal vector of Let [/] = — f- where f± is the one side limit of any quantity 
/ on r. For any test function (p £ Cq"(M^"^^) one may compute that. 



where V denotes the spatial gradient and the matrix R := diag(l,l,— 1) performs a 
reflection across the xt-plane. By taking cp identically zero on F, we obtain that, if u 
has a jump discontinuity along F, then 

(CI) Nf = N2 or equivalently F is parallel to x = y or x = —y which are character- 
istics of Oxy 

Therefore, the surface takes the form of 



2. Weak solutions and Compatibility Conditions 




idtu + OxyU + |npu = 0, 





T = {{t, xo{t) + ais, yo{t) + a2S \ t, s 



G R} and Nq = No{t) 
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where qi^2 = =tl- For such smooth F, integrating by parts on the first integrand leads 
to 

(2.3) j (pRN ■ {i[u],2[Vu])dS = 

and thus we have 

(C2) iMiVo + [nJiVi - [uy]N2 = 0. 
Condition (C2) is equivalent to that [ux], [uy], i[u] are colinear in the complex plane. 
In fact, if u is piecewise and [u] = 0, then, on the one hand, (C2) implies that 
the vector (A'^i, —N2) lies in the kernel of the 2x2 matrix [Vu]. On the other hand, 
if [Vn] 7^ 0, then its kernel is given exactly by the tangent space of F, i.e. the one 
dimensional subspace spanned by {N2, —Ni), and thus A''^ = N2 holds. This means 
that, even though (CI) follows from the assumption [u] 7^ 0, it still follows from (C2) if 
[u] = and [Vit] 7^ 0. Therefore, (CI) and (C2) should be imposed on weak solutions 
as long as u or Vn have a jump discontinuity along T. Condition (C2) can also be 
viewed as a necessary condition for the singularity to propagate as a surface in R^+^. 
If (C2) is not satisfied by some given initial data, then the singularity will not maintain 
a simple geometry for t > 0. 

Notice that, according to (CI), the jump discontinuity can occur along lines in the 
directions of (1, 1) or (1, —1). To study the local behavior of the weak solutions at the 
intersection of its lines discontinuity, we consider the case where 

= {{x - xo{t)f = (y - yoit)f}, 

i.e. a cross in the xy-plane with a moving center c{t) = (xo(t), yo(0)- We will use 
u^^ to denote the value of u in the cone either in the horizontal or vertical directions, 
opening in the positive or negative directions. In particular, u^'^^'^ denote the limiting 
values of u^^ on F in each quadrant. We use A'^ to denote the unit outward normal 
of the horizontal cones (x — XQ{t))'^ > {y — yo{t))'^. In this case, instead of ()2.3p . 
integrating by parts on the first integrand in (12. 2p implies 

4 j 0(c(t)){4(c(t)) - ulicit)) + ulicit)) - ulicit)) + uUc{t)) - ulicit)) 

(2.4) + ui{c{t)) - ui{cit))}dt - J (pRN ■ {i[u],2[Vu])dS = 0. 

Therefore, in addition to (C2), we also obtain 

(C3) ul-ul + ul-ul + ul-ul + ul-u'^ = along c{t), t G M. 

Remark 2.1. When a solution has singularity along a surface, the compatibility con- 
ditions ( C1-C3) above are only necessary. In fact, one should be careful when talking 
about the propagation of the singularity as the system does not have finite propaga- 
tion speed, which can be seen from the formula of the fundamental solution if the free 
equation |18] . For example, if the initial data has some simple jump discontinuity, this 
singularity might become unpredictable for t > except for highly symmetric cases. In 
what follows F will describe the cross X in defined by { |x| = \y\ }. We will see 
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in Section 4 that the simplest example of data satisfying the compatibility conditions 
described above are those of the form {f,g, — /, —g) given in each of the four cones of 
R'^\X. 

3. Hyperbolic coordinates and radial solutions 

Notice that the regular Euchdean rotation does not keep invariant, and hence 
the usual polar coordinates may not be the most appropriate ones for this system. 
Instead, we introduce an alternative: we foliate \ X by hyperbolas and consider the 
following change of variables: 



(3.1) 



X := a cosh a y := osinha := — when |x| > \y\ 

X := asinha y := a cosh a := — when \x\ < \y\ 



where o, a G M, a 7^ 0. 

Under this change of variables we have that the area form dxdy = \a\dada. Moreover 
for any a G M \ {0} we have: 

X y 1 
Va = (— , — ) and 00=- for |x| > 

a a a 

X y 1 

Va = ( , — ) and Da = — for |x| < |y|. 

a a a 

For Va we need to distinguish whether a > or a < 0. We have: 



Va = 


(-- 


X 


) for 


x\ 


> 


\y\ 


and 


a 


> 


Va = 




X 


) for 1 


x\ 


> 


\y\ 


and 


a 


< 


Va = 




X 


) for 


x\ 


< 


\y\ 


and 


a 


> 


Va = 




X 


) for 


x\ 


< 


\y\ 


and 


a 


< 


□a = 





in all cases. 















We moreover note that 

— tty = 1 and — = ^ for |x| > \y\, 

— tty = —1 and — = — for |x| < \y\. 

By changing coordinates (x, y) — > (a, a) in M2\X the (3HNLS) equation transforms 
itself into the system [| 



2 Note that in the first quadrant a > and a > 0; in the second quadrant a < and a < 0; in the 
third a < and a > and in the fourth quadrant a > and a < 
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(3.2) 



jidtu + daaU + ldaU- ^daaU+\u\'^u = for > |y| (F) 
[idtu - daaU - ^daU + ^daaU+\u\'^u = for < |y| (D) 

It is clear that (j3.2p is invariant under the translation in a, or equivalently (3HNLS) 
is invariant under the rotation in the hyperbolic coordinates. Thus it motivates us 
to study radial or equivariant solutions. We will focus on radial solutions u{t, a, a) = 
u{t,a). (See Remark 13.91 for equivariant solutions.) In the rest of this note, we will 
consider radially symmetric solutions assuming u^^ is a priori smooth for a G [0, oo) 
except for some jump discontinuity across X. Note in what follows these facts will be 
used without any further mention. In other words we will be working with 



(3.3) 



{idtu + daaU + ^daU + \u\'^u = for |2;| > |y| (F) 
idfU — daaU — \daU + \u\'^u = for l^l < \y\ (D) 

The compatibility condition (C3) implies that at the origin (0, 0) G 
(3.4) ^+(o) + n-(0)=n+(0)+<(0). 

Some obvious solutions of this type are {u^ = —u^,u^ = = 0) or {u^ = = 
0, u~ = — n+) where or n+ solves the above (F) or (D) respectively. 

We are interested in the existence (construction) of radial standing wave solutions 
u{t,a) = e~^^^v{a) such that as a — )• ±oo v — >• v±oo with v±oa a constant. 

Remark 3.1. Ifu(x,y) is and (hyperbolically) radially symmetric, then the radial 
symmetry implies Du{0,0) = and thus Va{0) = 0. We assume the latter from now 
on. 

From (13.31) we then have that v satisfies: 



(3.5) 



l^v + Vaa + + \v\'^v = for |x| > |y| (F) 

fJ'V - Vaa - Iva + \v\'^V = for |x| < (D) 



Remark 3.2. We callvd the solution to (D) andvf the solution to (F) andv = {vd,Vf) 
the solution to ()3.5p . Note that a radial solution to ()3.5p is if and only if and 
Vf have equal values at a = 0; ie. t'rf(O) = Vf{0). 



Remark 3.3. Under the assumption Va{0) = 0, without any loss of generality we 
consider v real valued only. In fact, suppose v{0) = woe*^" with omegao G K and 
6q G M. Let v{a) = v{a)e~^^° , then v{0), Va{0) G M. Theorem \7.1\ implies the solutions 
v{a) exist and are real. 



Remark 3.4. In the analysis below we only consider a > but corresponding results 
hold for a < 0. 
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Remark 3.5. By considering another coordinate change x' = a cos a and y' = a sin a, 
(j3.5p (F)(D) simply become ^/S.xiy'V+ ^v+v"^ = 0. Solutions of this equation in L'^ have 
been studied extensively (see, for example, [U [21 [3l [261 [27] j. Based on those results, it 
is obviously impossible to find solutions of (F) and (D) with decay at infinity so that 
they match at a = 0. Instead, we consider all possible bounded solutions here through 
a more careful analysis. 



Theorem 3.6 (Main Theorem 1). weak radial standing wave solutions u{t,a,a) = 
e~*^*u(a) to (j2.ip . which are except along = (or equivalently a = 0), exist 
only if fi < 0. In fact, for any vq £ (0, there exists a unique such solution 

{uf,U(i) = (e~*'^*vj(a), e~*^*Wrf(a)) such that Vf and v^ are smooth on a £ [0,oo) and 

Vd{^) = f/(0) = vo and daVd{0) = = daVf{0). 

Moreover, for a — )• oo, 

Vf(a) = \/^u H j= cosha + aloga + /3) + 0(a~2) 

Vdici) = —F= cosHa + crloga + /3) + 0(a~2) 
for some 7, /3, o" G M and 7^ and r^ 7^ 0. 

For Vq G (— ^, 0), we only need to consider —v and apply the above theorem. 

Notice that equation (3HNLS) is invariant under the scaling transformation u — >■ 
u\{t, x, y) = \u{\^t, Ax, \y). Accordingly, in the hyperbolic coordinates, consider (j3.2p 
in the usual similarity coordinates ii(i,a,Q;) = t~^u{[og f~2a, a), we obtain 

(3.6) 

iidtu - f -u - ^adaU + daaU + \dau - -ijdaaU + \u\'^u = for |x| > \y\ (F) 
[idtu - ^u- ^adaU - daaU - \daU + -^daaU + \u\'^u = for |x| < \y\ (D) 

where now u = u{t,a,a). Next we look for radial standing waves in the self-similar 
variables of the form u = e~^^^v{a) which satisfies 



(3.7) 



^j.v - ^v - ^aVa + Vaa + \va + \v\'^v = ^ for |x| > |y| (F) 
^iv- \v - ^ava-vaa-\va + \v\'^v = for |x| < |y| (D) 



The remarks we made for the radial standing waves are valid for radial self-similar 
waves except here we can no longer assume v has only real values. 

Theorem 3.7 (Main Theorem 2). For any /i £ R, there exists f * > such that for any 
Vq G (— i;*, u^k), there exists a unique weak radial self-similar wave solution {uf,Ud) 
of (j2.ip which is except at x^ = y^ ( or equivalently a = j and 

Uf{t,0) = Udit,0) = ?Jot"^e"*^^°sl*l and dau{t,0) = 
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As ^ oo, there exist constants rj / 0, / 0, 9f, and 9^ such that 



2 - 
u^{t,a) = !:^e^(f^-Mi°gl*l) cos(% - ^ - 2/. log ^) + 0{%) 
a at yjt a-^ 

2 - 

u,(t, a) = !i^e^(-|^-A'i°g 1*1) cos(0, - ^ + 2/. log 4) + O(^). 
a or yjt a 



Remark 3.8. From the asymptotic form of the above self-similar waves, for any a 7^ 0, 
i.e. away from the cross \x\ = \y\, as t approaches the singular value (here t may be 
replaced by ib(t — to)), the complex values of the solutions oscillate rapidly in both the 
angle and the magnitude. In particular, their magnitudes oscillate with an envelope 
proportional to ^. At the cross \x\ = \y\, the magnitude of the solutions blows up like 
. Therefore they strongly display a kind of profile of the X shape. 

To prove these theorems, we first prove existence of solutions near a = and next we 
study their properties as a — )• +00. Equations (|3.5p (F)(D) and (|3.7p (F)(D) have been 
extensively studied. To establish local existence of solutions near a = 0, well known 
methods include for example, power series or variational methods. Here we present (in 
the Appendix) another method based on invariant manifold theory to study solutions 
near a regular singular point of a class of ODE system of which (13.5[) (F)(D) and 
()3.7p (F)(D) are special cases. The advantage of this method is that it does not require 
either much smoothness as in the power series approach or variational structures. More 
specifically, for any vq, let xi = v — vq and X2 = ava in these four ODEs. Each of them 
turns into a form for which we may apply Theorem l7.1[ For example, for (j3.5p (F). in the 
notation of Theorem 17. II we have that F(a, x, vq) = (x2, —^a?{xi +vq) — a^(xi +^0)^)"^ 

with the parameter vq and -F(0, 0,fo) = for all vq. Moreover Fa;(0,0, uq) = 

and thus Ej^ = {0}. Therefore, for any vq, there exists a unique solution v[a) for small 
a > so that v and aVa are Holder continuous in both a and v{{)) = vq. Moreover, the 
statement (4) in Theorem 17.11 implies X2 = aVa = O(a^) which implies Va{0) = 0. The 
other equations are handled similarly. 



Remark 3.9. The same Theorem 7.1 can also be applied to study equivariant standing 
waves or equivariant self-similar solutions where the equivariance means u{t, a, a) = 
e^'^e~^^^v{a), c, /u G M. In this case, for equivariant standing waves, (13. 2p implies 

( fJ-v -\- Vaa + lva-\- -\-\v\'^v = for \x\ > \y\ (F) 

1 fJ-V - Vaa - ^Va - -\- \v\'^v = for \x\ < \y\ (D) 



For c 7^ 0, Theorem \ 7.1\ immediately implies that the only continuous solution on 
[0,ao), oq > 0, of either of the above has to be trivial. Therefore, there are no nontrivial 
continuous equivariant standing waves. The same argument applies to the self-similar 
case as well. 
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To finish the proof of the theorems, we will study the asymptotic properties of these 
systems as a — t- oo in the next two sections. 



4. Radial Standing Wave Solutions 

Note that by considering oi = |//|2a and v = |^|2fi, we only need to consider 
IJ: = ±1,0 for (j3.5p fF)(D). To study the asymptotic properties of the solutions as 
a — >• +00 we proceed according to the sign of fi. 

Case ^ > 0. First we show that all solutions vj to (|3.5p (F) go to zero like cos{'ja+ 
a log a + (3) as a — )• oo. Indeed, for nontrivial v satisfying Vaa + + fJ.v + \v\'^v = 
and a > define v =: -^w. Then w satisfies: 



(4.1) 



/I ,1 1 1 1 , ,2 

+ fJ-j—^W + —^Waa H l=\w\ W = 0. 

4a^ \/a \ a a \ a 



We want to first show that w is bounded. Multiply (|4.ip by i/a to get 



(4.2) 



/I N 1, ,2 

Waa + K^r^ + [i')w + - li; = 0. 

4a^ a 



Multiplying (14. 2p by ^^;a and rewriting the resulting expression using 'perfect deriva- 
tives' we get 



(4.3) 



da 



2^) + ^- +2- +4^- 



1 

'4^ 



w 



1 

4^ 



Define 
(4.4) 

Then we have 



M 2 



1 



8a2 



+ -w^ + — ii;^ > 
2 4a ~" 



(4.5) 



a.i.(a) = -^u;^-^^.^<0 



Hence i?(a) is decreasing, which implies w and are bounded and Eq := lima_!.oo -£-(0) 
exists. The boundedness of w and (14. Sp again imply 

> > -^E 

for some C > independent of a. Therefore £"0 > and E{(i) — Eq = O(^). Moreover, 
let us change the variables 



w = r cos a, Wa = r sm ( 
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Recall = 1, we have — 2Eo = O(^) and thus 



(4.6) r = V2^ + 0(-) 

a 

as a — )• oo. One may compute 

1 1 s 

Ta = — —^r COS sin 9 r cos 6* sin 9 

4a^ a 

= - 1 - i cos^ 9 - -r^ cos^ 9 
4a^ a 

3i?o I -E^o / 3 4 -Eo - 4 

= -1 \ (--cos 0)-\ cos 0. 

8a a 8 a 

Integrating the 9 equation by parts and using 6a = —1 + O(^), one may show, 

(4.7) = 0^ + a + ^\oga + Oi-). 



We next consider solutions Vd to (|3.5p (D). Multiplying (D) by av and integrating 
on [0, a], we obtain 

{v\ = l f ~a{{vaf + v^ + v^) d~a 

,,2 







which implies v > O(loga) as a — )• cxd unless v = 0. Therefore there are no nontrivial 
bounded solutions Vd to (j3.5p (D). 

All in all we conclude then that for the case // > there are no bounded solutions 
other than (w/, 0, —Vf, 0), which can not be in unless trivial. 



Case n = 0. For (|3.5p (F), let =: a one may compute 

1 1 _2 ^ 

Waa + ir'^a + TT^'^ + ^ = 0- 

3a Oa"' 

2 

Multiplying it by a^Wa, we have 

9a(-a3(u'a)^ + -w* + — a 3u;^) = -— a 3^^;^<0. 
z 4 18 2 / 

Therefore w = 0(1) and Wa = 0(a~3) as a — 00, which implies 

2 1 
t'a = 0(a~3), v = 0{a^3), as a — )• 00. 



In fact, for different values of f (0), the solution v{a) only differs by a scaling change. 

Finally, the same argument as in the case of /u > shows that the only bounded 
solution Vd to ()3.5p (D) is the trivial one. 

We thus also conclude in the case /i = that there are no bounded solutions other 
than (f/,0, — fj,0), which can not be in unless is trivial. 
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Case fi < = — !)• From the uniqueness given in (4) of Theorem 17.11 of solutions 
which are Holder continuous at 0, clearly Vd,Vf = 0, ±1 if ?j^(0), ?;/(0) = 0, ±1. 

We first consider solutions to (j3.5p (D) and claim that if |frf(0)| > 1, then > 
1 and — )• oo as a — )• oo. Indeed, without any loss of generality assume Vd{0) > 1 
and Vd{a) = 1 + z{a). One may compute 

(4.8) Zaa + -Za - 2z - 3z^ - Z^ = 0. 

a 

Multiplying it by Za, we have 

2 4 a 

which implies 

> -{z{Of + z{Of + lz(0)^) > (^(z.(g))2 - {z{af + z{af + iz(a)^)). 

Therefore, z{a) / 0, for all a > which yields z{a) > 0. Multiplying (j4.8p by —az 
implies 



da{z') = - 

a 



i-a 

/ a{{zaf + 2z^ + 2,z^ + z'^) da 
Jo 



which implies > clog a. 



Thus we assume Vd{0) G (—1, 1) and claim Vd ^ like cos(7a + crlogo + /3) as 
a — )• +00. Indeed, multiplying (j3.5p (D) by Va we obtain 

2 2 4 a 

Since — ^v'^ is strictly increasing in for G [0, 1], we have 

(4.9) v{0) G (-1,1) ^ v{af <v{Of <1, a > 0. 

We will follow closely the procedure which we used to handle solutions to (|3.5p (F) in 
the case of u > 0. Define v =: A=w. Then w satisfies: 

(4.10) Waa + ( A + --W^ = 

4a^ a 

Multiplying (j4.10p by Wa and rewriting the resulting expression as 'perfect derivatives' 
we get 



-(■w„) -\ H — w w 

2^ 8o2 2 4a 



4a3 4a2 



(4.11) da 
Define now 

(4.12) Eia) := i(u;^) + + ^w' - j-^ > ^(u;^) + (i - M2)^)^^ 



where we used \a 2w\ = \v\ < \v{0)\. From (|4.10p we have that 
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Using |a~2t[;| = \v\ < \v{0)\ and (14. 9p again, we obtain from ()4.13p 

a ' 2-|i;(0)|^ 



(4.14) daE{a)<-E, c= J ^1,2 < 1- 



The differential inequahty (|4.14p imphes E < 0{a'^) and thus w < 0(a2 ) as a — >■ oo. 
Now from ()4.12p and (|iT9|) we obtain that 

(4.15) E{a) > Kow"^ where ^ > := ^ - ^^^^^^ > 0. 
Using (|4.13p and (|4.15p we conclude that 

(4.16) daE{a) < 0(a"-2)w;2 < o{a^-'^)E. 

Hence E = 0(1), which gives w = 0{1) and Wa = 0{1) from (|4.12p . The oscillatory 
asymptotic form of v then fohows from the exactly same argument used in the case 
/i > 0. 

Next, we consider solutions Vf to (13. 5p (F). We will show that for f/(0) G (0,2), 
solutions Vf to (j3.5p fF) behave like 1 + -^cos(7a + uloga + /3) as a — )• +oo. By 
the odd symmetry of the equation, the case when f/(0) G (—2,0) is the same. First, 
multiplying p.5p (F) by Va, we have 

which, along with Va{0) = implies that 

(4.17) \v{af -l\ <\v{Of 

In particular, (|4.17p yields v{a) > 0. Again, let v =: 1 + a^^w; then w satistfies 

(4.18) Waa + {2 + ^)w + 3a-iw^ + -w'^ = 0. 

4a^ a 

Multiplying (j4.18p by Wa, we have 

(4.19) daE = -^a~^w^ - ^a~^w^ - ^a'^W^ = -^a"^ii;^(a"^ + 2a~^w + a~^w'^) 

where 
(4.20) 

In the last two inequalities, we used the fact v{a) > 0. For the same reason, we can 
estimate the right hand side of of ()4.19p as 

2 

daE <^(l-v'^) <-E, c = \v(Of - l\ < I. 
4a a 

It implies E = 0{a^) and thus w = 0(a2) as a — )• oo. Substituting this into (|4.19p we 
obtain 

daE = 0{a-^)E 
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and thus E = 0(1) and in fact E^EQ>Oasa^oo unless w = 0. The oscillatory 
asymptotic form of w can be obtained much as in the previous cases. 

In fact, by a more careful analysis, one may show the behavior of at a — )• oo 
is always one of oscillatory convergence to ±1 except maybe for a sequence of data 
Vf{0) = Vj, \vj\ > 2 which make Vf converge to exponentially. Since we seek bounded 
continuous solutions to (|3.5p and we have already showed that solutions Vd to (j3.5p (D) 
are unbounded if |fd(0)| > 1, the above result is sufficient for us. 

Prom the above we conclude then that for the case // < there exist smooth and 
Vf so that Vd{G) = Vf{0) G (—1, 1), daVd{0) = = daVf{0), Vf{oo) = ±1 and Vd{oo) = 
where the convergence is like a~ 2 cos (7a + (3 log a + a + 0{^)). 



5. Radial self-similar waves 

To analyze the system (|3.7|) for self-similar solutions, we first change the independent 
variable to s = ^ and then change the unknown to w = s'ie^i'^v. Consequently, (j3.7p 
becomes 



(5.1) 



Wss + ije + ^ + ^,)w + ^\w\^w = for |x| > lyl (F) 
Wss + ijQ + ^ - - ^\w\'^w = for l^l < (D) 



Without loss of generality (see Remark |3. II and [3^ . we may assume w is real- valued. 
It is clear that, as long as we obtain the boundedness of the solutions to the ODEs, 
their oscillatory asymptotic forms follow from the exactly same change of variables as 
in the previous section w = rcosO and Ws = \rsm9 and similar arguments. So we 
focus only on the boundedness. 

Multiplying (|5.ip by Wg and letting 

(5.2) E{s) = -wl + (— + ^ ± ^)w'^ ± -^w^ 

we obtain 

which we would like to control in terms of (|5.2p . 

Focusing equation (F). In this case, for s > sq := 16(|/i| -|- 1) we have E > and 

^) j?^ \ ^ 1^1 2 ^ MM ^ 



Therefore E{s) — t- i?o > as s ^ cxd which implies that w and Ws are bounded as 
desired. 
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Defocusing equation (D). In this case, we follow a standard bootstrap argument. 
For those solution w{s) such that |tf(so)| < 1; let 

si := sup{s > So I |w(s')l ^ l^'s' £ [soj-s]}- 
For s G [so, si], we have from (15. 2p 

(5.4) E{s) > ^w', + ^w' 

and from (15.311 and (15.411 we obtain that 



Therefore, if |u)(so)| and |'Ws(so)| are sufficiently small -which by the continuity in 
Theorem 17.11 is guaranteed by taking w{0) sufficiently small- we have that si = oo. It 
follows w{s) and Ws{s) are bounded. 



6. Conclusions & Future Directions 

In conclusion, in the present work, we have considered the two-dimensional hy- 
perbolic NLS equation. We have illustrated the relevance of the use of the so-called 
hyperbolic variables in the context of this PDF model (relevant to optical media com- 
bining normal dispersion -typically in the time variable- to anomalous one -typically 
in the continuous spatial variables). We have discussed the context of weak solutions 
within the model and how relevant compatibility conditions naturally arise along the 
characteristic lines of the hyperbolic operator. We have subsequently elucidated the 
prototypical weak standing wave solutions that the model supports establishing that 
they only exist for negative values of the propagation constant fj, and characterizing 
their weak modulated power law decay on the basis of energy-type methods (in the 
radial hyperbolic variable). Also, similar techniques but in appropriately rescaled vari- 
ables have been used to discuss the existence of radial, self-similar wave solutions to 
the PDF for all fieM.. 

It would be particularly interesting to revisit the numerical computations of earlier 
works [11^ [131 fT4l \TE[ [TT] in the context of weak solutions and to examine whether 
for appropriately crafted initial conditions the relevant direct numerical computations 
support the decay rates analytically obtained herein. Preliminary computations with 
radial initial data such as u{x, y, 0) = tanh(a) for |x| > \y\ (and zero otherwise) seem to 
support the radial evolution of the solution profile and the development of oscillations 
around the unit asymptotic state, as suggested by our Main Theorem (see Fig. [T]l. In 
particular, this initial profile provides a vanishing solution in the (D) region, as well 
as a radial solution asymptoting to unity in the (F) region. The observation that can 
be made based on these preliminary results is that in the latter region the data indeed 
remains radial in nature (since the hyperbolae forming in the middle panel are equi-a 
lines) yet it also develops an oscillatory dynamics associated with the decay to the 
stationary state, which is reminiscent of the oscillatory convergence suggested by our 



HYPERBOLIC NLS 



15 



Main Theorem. It should also be noted that the numerical method used here was a 
4th order Runge-Kutta for marching the system in time, combined with a centered- 
difference in space scheme in a sufficiently fine spatial grid. However, detailed numerical 
investigations are relevant to establish the relevant decay; these will be deferred to a 
future publication. Prom the point of view of analytical considerations, it would be 
relevant to generalize the results obtained herein to the more experimentally tractable 
directions of 3-dimensional media (with two focusing and one defocusing direction) as 
in |13 1 114 1 [TT t 112). In that setting, it would be relevant to connect generalizations of the 
present results to the rapidly evolving literature on the nonlinear X- waves of the above 
works. Another highly promising direction is that of the quasi-discrete waveguide 
array setting of equation (jl.ip in 2 + 1-dimensions or of |22j in 3 + 1-dimensions. 
Notice, however, that even the setting of equation (3HNLS) considered herein may be 
of relevance to physical applications of light propagation within planar waveguides in 
glass membrane fibers, where already some interesting phenomena such as linear and 
nonlinear guidance with ultralow optical attenuation have been observed; see e.g. [23]. 

7. Appendix: An ODE blow-up theorem near a regular singular point 
Consider 

(7.1) x' = ^F{t,x), t>0,x£W, F e C7'=(M"+\R"). 

Here t = is a so called regular singular point of the system and we are interested 
in solutions continuous or smooth at t = 0. A subcategory of this system is a scalar 
valued equation 

which include (j3.5p (F)(D) and (|3.7p (Fl(D) as special cases. In fact, by setting Xj = 
t^~^(y^^~^^ — Cj) for j = 1, . . . ,n and arbitrary constants ci, . . . , Cn, the above equation 
turns into the form of (|7.ip . 

Let B(X,r) denote the ball centered at and with the radius r in a Banach space 

X. 

Theorem 7.1. Assume k >1. 

(1) There exists a solution x{t) to (|7.ip which is continuous on t & [0,e] for some 
e > and x(0) = xq if and only if F(0, xq) = 0; 

(2) Suppose F{0,xq) = 0. Let be the generalized eigenspace of Rj. {0,xq) cor- 
responding to all eigenvalues whose real parts are greater than 0. Then there 
exist e > and a C'^ mapping h : B{E^,e) — B{{E^)-'- , e) so that: x{t) is 
a solution of (j7.ip on (0,e] with x{e) G B{E^,e) + B{{E^)-^ ,e) and Holder 
continuous for t G [0,e] and x(0) = xq if and only if x{e) G graph{h). 

(3) The above solutions are of order 0{t^) where /3 G (0, /3o) and 

/3o = min{l, inf{Re A \ X is an eigenvalue o/Fx(0,xo) and ReA > 0}}. 
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Figure 1. Evolution of Eq. (|2.ip for 7 = 1 (a factor of 1/2 is also 
used in front of the as is often done in the physics literature). 
The top panel shows the initial condition and the middle panel shows 
the result of the evolution at t = 1. Both panels show contour plots 
of I up. Notice the hyperbolic equi-a contours forming in the middle 
panel. The bottom panel shows a radial cross-section (at y = 0) of 
the dynamics, illustrating the oscillations around the asymptotic state 
that spontaneously develop in the dynamics (as suggested by our Main 
Theorem). The dashed line shows the initial density profile, the dash- 
dotted one shows its evolution at t = 0.5 and the solid one its evolution 
at t = 1. 
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(4) Ifk>2 and Ej^ = {0}, there exists only a unique solution x[t) on (0,e] so that 
it is Holder continuous in t on [0, e] and x(0) = xq. Moreover, if Ft{0,XQ) = 0, 
x(t) = 0{t^) ast^O. 

(5) If F = F{t,x,a) with a external parameter a, F ^ , and F{0,xo, •) = 0, the 
mapping h in item (2) is also C^~^ in a. 

Proof. To blow up the system to analyze solutions near t = 0, we create an auxiliary 
independent variable r so that fp = t, use X to denote the differentiation in r, and 
rewrite (|7.ip as an autonomous system 

(7.2) I i - i or i = Fix) 

\x = F{t,x) 

where x = (t,x)^. A solution of ()7.ip is continuous at t = is equivalent to that its 
corresponding solution of (|7.2|) converges as r — )• — oo. On the one hand, it is clear 
that xo = limT-_!._oo 2;(r) only if xq = (0, xq)^ is a fixed point of (|7.2p . On the other 
hand, suppose -F(0, xq) = 0, we notice 1 is an eigenvalue of DF{xo) with a generalized 
eigenvector T = (l,w)^, for some w G M", and -E+ = £^+©span{T} is the generalized 
eigenspace of DF{xq) corresponding to all eigenvalues whose real parts are greater 
than 0. From the standard unstable manifold theorem (c.f for example [4l[5l[6]), there 
exists a manifold C R"+-^ in a neighborhood of xq so that (a) xq G W^] (b) 
TxqW^ = £'+; (c) it is locally invariant under the flow of (|7.2p ; and (d) x G if 
and only the solution x(t) of (I7.2p with x(0) = x stay close to xq for all r G (— oo,0] 
and x(r) — >• xq exponentially as r — >• — oo. Changing the variable back from r to t, 
the exponential rate 0{e^'^) as r — )• — oo becomes 0{V). Therefore, for equation (|7.ip . 
solutions on W'^\{t = 0} exactly correspond to solutions which are Holder in t and 
converge to xq as t — t- 0. This proves statement (2). Statement (3) is obvious when 
one notices that t = ce^ . Statement (4) is a consequence of the tangency of to 
which is equal to MT = R(l, 0)^ in this case, at xq- Statement (5) follows directly from 
the the standard result that unstable manifolds are smooth in the external parameters 
El. □ 
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